High-fidelity readout of trapped-ion qubits 
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We demonstrate single-shot qubit readout with fidelity sufficient for fault-tolerant quantum com- 
putation, for two types of qubit stored in single trapped calcium ions. For an optical qubit stored 
in the (4Si/2, 3D 5 / 2 ) levels of 40 Ca + we achieve 99.991(1)% average readout fidelity in one million 
trials, using time-resolved photon counting. An adaptive measurement technique allows 99.99% 
fidelity to be reached in 145 /is average detection time. For a hyperfine qubit stored in the long- 
lived 4Si /2 {F = 3, F = 4) sub-levels of 43 Ca+ we propose and implement a simple and robust 
optical pumping scheme to transfer the hyperfine qubit to the optical qubit, capable of a theoretical 
fidelity 99.95% in 10 fis. Experimentally we achieve 99.77(3)% net readout fidelity, inferring at least 
99.87(4)% fidelity for the transfer operation. 



A quantum computer (QC) requires qubits which can 
be prepared and measured accurately, and made to in- 
teract through high-quality quantum logic gates. Suc- 
cessful fault-tolerant operation requires certain minimum 
thresholds for the fidelity of state preparation, logic gates 
and state measurement. Qubit readout is vital, not 
only for the final output of the QC, but for the error- 
correction essential for its operation. There is a trade-off 
between the error rate permitted and the number of extra 
qubits required to achieve fault-tolerant quantum error- 
correction (QEC), but typical studies require errors to be 
below ~ 10~ 3 for realistic implementations [H, Q. A cor- 
rection network typically uses more gates than individual 
qubit readouts, so readout error is less critical than gate 
error [2]. On the other hand, precise readout can be 
used to compensate for gate error. It is also crucial in a 
measurement-based QC [3]. 

Trapped-ion QC approaches generally use qubits based 
either on a hyperfine transition within the ground atomic 
level or on an optical transition between ground 

and metastable states |{|. State preparation is im- 
plemented by optical pumping. State measurement is 
achieved by repeatedly exciting a cycling transition which 
involves one of the qubit states but not the other, and 
measuring whether or not the ion fluoresces [Toj ] . The 
accuracy with which the measurement can be performed 
depends on the rate at which fluorescence photons can be 
detected for the "bright" qubit state, compared with the 
rate at which the "dark" state gets pumped to the bright 
manifold (e.g. by off- resonant excitation), or vice versa. 
The fact that many fluorescence photons can typically be 
emitted by the ion before the undesired pumping process 
occurs allows for high-fidelity single-shot measurement 
even though the absolute photon collection efficiency can 
be poor (typically ~0.1%). Measurement fidelities of 
98%-99% have been reported 0, [1,0,11. In recent work, 
an optical qubit was measured indirectly by repetitive 
measurement on an ancilla qubit held in the same ion 
trap, yielding a fidelity of 99.94% in — 12ms [ill ]. 

In this Letter, we report direct high-fidelity measure- 
ment of firstly an optical qubit stored in the (4Si /2 , 3D 5/2 ) 



levels [18[ of a single 40 Ca + ion (fig.[U right) and secondly 
a hyperfine qubit stored in the S 1/2 (F = 3, F = 4) sub- 
levels of a 43 Ca + ion (fig. [3j right), where we first map 
the hyperfine qubit to the 43 Ca + optical qubit. Readout 
is achieved by driving the (Si /2 <-> Pi /2 <-> D 3/2 ) mani- 
fold and detecting the fluorescence from the P 1/2 — > S 1/2 
decay: fluorescence indicates the qubit was initially in 
the "bright" Si /2 state, absence of fluorescence that the 
qubit was in the "dark" metastable D 5/2 state (lifetime 
t = 1168(7) ms [l2|). We discuss state inference by sim- 
ple photon-count thresholding and by maximum likeli- 
hood methods using time-resolved detection. The latter 
are applicable to a wide range of physical qubits [l]| • 

For the Ca + measurements, a single ion is held in a 
Paul trap [T3|, Doppler-cooled on the Si /2 <-> P 1/2 transi- 
tion by a 397 nm laser and repumped on D 3/2 <-> P 1/2 by 
an 866 nm laser. A photomultiplier (PMT) detects fluo- 
rescence from the ion with a net efficiency 0.19(2)%. The 
ratio of ion fluorescence to background scattered laser 
light has a maximum value « 690 at low 397 nm intensity, 
but optimum readout fidelity is achieved at higher inten- 
sity because the increased fluorescence can be detected 
more rapidly compared with r. Mean photon count rates 
for fluorescence (i?s) and background (Rd) in this exper- 
iment were Rb = 55800 s -1 and Rd = 442 s" 1 (including 
the PMT dark count rate of 8.2s -1 at 20°C). 

The experimental sequence consists of preparing and 
measuring each qubit state in turn, and comparing the 
measurement outcome with the known preparation. To 
prepare the D 5/2 "shelf" state a 393 nm beam is switched 
on to drive S 1/2 <-> P 3/2 for 1 ms; spontaneous decay 
populates D 5/2 (the 866 nm beam is left on to empty 
D 3/2 ). The rate of transfer to the shelf was measured to 
be [12(2) /us] -1 . The 393 nm beam is extinguished and, 
20(2) /Lts later, we start to collect PMT counts. Photons 
are counted for a bin time of 2 ms, which is divided into 
200 sub-bins of duration t s = 10 /is. Two alternately- 
gated 10 MHz counters are used, so that there is negli- 
gible (< 50 ns) dead-time between sub-bins. The counts 
rii from each sub-bin i are available to the control com- 
puter ~ 5 /is after the end of the sub-bin, should they 
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FIG. 1: 40 Ca + photon-count histograms with t& = 420 ^s, for 
bright S1/2 state (o) and dark D 5 / 2 state (•) preparations. 
Above-threshold events in the dark histogram give en ; below- 
threshold events in the bright histogram give e_e. Insets show 
time-resolved counts for two individual 17-photon events, one 
from each histogram, with likelihood ratios Pb/pd- A PMT 
dark count histogram is also shown (a), normalized to the 
same area. Dashed curves are Poisson distributions with the 
same means as the histograms. The solid curve is the expected 
D 5 / 2 distribution taking into account spontaneous decay. 



be required for real-time processing, and are recorded in 
time order. Next, the ion is prepared in S 1/2 by driving 



D5/2 P3/2 for ~3ms with an 854 nm laser beam, after 
which fluorescence is collected for a second 2 ms bin. 

This sequence was repeated to give 2 20 w 10 6 trials, 
in half of which the ion was prepared in the dark D 5/2 
state and in half of which it was prepared in the bright 
S 1/2 state. We define the average readout error to be 
6 = \{ € b + £d)j where es is the fraction of experiments 
in which an ion prepared in the bright state was detected 
to be dark, and similarly for cd- In the case of the dark 
state experiments, there are two small preparation errors 
arising from (i) the finite shelving rate of the 393 nm 
laser and (ii) the probability of spontaneous decay during 
the 20 [is delay between state preparation and the start 
of the detection period. These give a contribution to 
e D of (12 (is + 20^s)/r = 0.28(2) x 10~ 4 . Known state 
preparation errors for the bright state are negligible, so 
the net contribution to e is 0.14(1) x 10~ 4 . Values of e 
given below are after subtraction of this quantity. 

Histograms of the number of photons n = Xh=i n i 
collected for the bright and dark state preparations are 
shown in fig. [TJ Here the first N = 42 sub-bins of the 
detection period were summed for a total detection time 
of tb = Nt a — 420 /its; this choice of tb and a threshold at 
n c = 5^ counts optimize the discrimination between the 
bright and dark histograms. The error e as a function of 
tb is shown in fig. (SJ for N = 1 . . . 100. As tb increases, 
e first drops rapidly due to decreasing overlap between 
the two histograms, to a minimum e = 1.8(1) x 10~ 4 at 
tb = 420 /is, then rises slowly because of the increasing 
probability of decay from D 5/2 during tb- 

At the optimum tb, the majority of the error is due 



to the above-threshold events in the dark histogram: 
erj = 3.2(3) x 10~ 4 . These mostly arise due to spon- 
taneous decay during tb, but we also observe above- 
threshold events from non-Poissonian PMT dark counts. 
A PMT dark count distribution is also shown in fig. [1] 
the highly non-Poissonian tail in the histogram probably 
arises as a result of cosmic rays 14j . The slow time-scale 
of luminescence in the PMT envelope excited by the cos- 
mic rays makes it difficult to eliminate these events en- 
tirely. We estimate that, for the threshold method, cos- 
mic ray events account for ~ 20% of erj. 

The threshold method does not make use of the arrival- 
time information of the photons. Using this information, 
we may hope to identify some of the events where decay 
from D 5/2 occurs during tb. For example, an ion decaying 
near the end of tb may give > n c counts but the fact that 
these occur in a burst at the end of the bin would suggest 
that this is more likely a decay event than a bright ion 
(compare fig. [T] insets) . The use of time- resolved mea- 
surement in the context of qubit readout was suggested 
in [1| and modelled theoretically in 1^, 15|; we follow a 
similar maximum likelihood treatment to Langer 15 1. 

We calculate the likelihood ps that a given set of sub- 
bins {rii} could have been generated by a bright ion, and 
compare this with the likelihood po that {n^} arose from 
an ion that was dark at the start of the detection period. 
If pb > Pd we infer that the ion was bright and vice versa. 
Pb = P({ni}|bright) is given simply by the product of 
probabilities ps — H^ =1 B(rii) where B(rii) is the proba- 
bility of observing ri{ counts from a bright ion (in the ideal 
case B(jii) is a Poisson distribution with mean count per 
sub-bin i?s^ s ). The calculation of pn = P({rii}|dark) is 
more involved, because we must sum over the possibil- 
ities that the ion decayed after or during the detection 
time tb = Nt s . To a good approximation 
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where (1—tb/r) is the probability the ion decays after the 
last sub-bin, (t s /r) the probability it decays during any 
particular sub-bin j, and D(rii) the background count 
distribution with mean Rot s - Using a recursion relation 
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N j-1 
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PD = (1- t b /r)M N + (t s /r)S N 

M = 1, M k = M fe _iZ>(n fc ) 

So = 0, S k = (S k -i + M fc _i)B(n fc ) 



(2) 



reduces the calculation of po from 0(A^ 2 ) to O(N) op- 
erations, making real-time readout much faster. In (|TJ) 
we made the simplifying approximations that t s ,tb <C t, 
and that the mean count rate changes from Rn to Rb at 
the start of the sub-bin j in which the ion decays; simi- 
lar expressions are found if these approximations are not 
made, with negligible effect on the results. 

The results of applying the maximum likelihood 
method to the same data set are also shown in fig. [21 
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for N = 1 . . . 100. For B(m) and D(m), Poisson distri- 
butions were convolved with the PMT dark count dis- 
tribution • We see that for > 200 p,s the maximum 
likelihood method gives a lower error than the photon- 
count thresholding method, tending to an asymptotic 
value e = 0.87(11) x 10~ 4 (with e D = 1.5(2) x 10~ 4 ), 
or a fidelity 99.9913(11)%. Furthermore, the maximum 
likelihood method does not require a precise choice of pa- 
rameters (n c ,tb) which must be determined from control 
data; we only need to know the independently measured 
rates Rb and Rd, and to choose a sufficiently long if,. 
The observed error is compared with that from a simu- 
lation of 10 9 trials using ideal Poisson statistics in fig. [2] 
although at short tb the experimental e is greater due to 
super-Poissonian noise in the data, at long if, it converges 
to the simulation's asymptote at e = 0.89 x 10~ 4 . We 
conclude that, at optimum if,, the maximum likelihood 
method is less sensitive to experimental noise (e.g. cosmic 
ray events, drift in Rb) than the threshold method. 

We now discuss faster detection by using an adaptive 
version of the maximum likelihood technique 1 1] . In the 
preceding discussion, the detection bin length tb = Nt s 
was fixed and we evaluated the likelihoods pb,Pd at the 
end of the bin, deciding that the ion was dark if pp> > pb- 
However, the absolute values of pb and pd also contain 
useful information. The estimated error probability eo 
that we have incorrectly deduced the ion to be dark when 
Pd > Pb is given by Bayes' theorem [l6j |: 



en 



1 - P(dark|{n. i }) = 1 



P({nJ | dark) 

P({m}) 



Pb 



Pb +Pd 



and similarly for es in the case where ps > Pd- By 
evaluating pu and ps at the end of each sub-bin, we can 
terminate the detection bin at t a < tb when the estimated 
error probability e# or eu falls below some chosen cut-off 
threshold e c . We also impose a cut-off time t c < tb in case 
the error threshold is not reached; t c is thus the worst 
case readout time. The result is that, for a given desired 
error level, the average detection bin length t a is shorter 
than when the bin length is fixed. If e c is not reached, 
es or ejj still quantifies confidence in the measurement 
outcome, useful for QEC. Using the recursion relations 
(|2"|) we can evaluate pb,Pd in < 1 /is, while each sub-bin 
is in progress, so there is negligible time overhead. 

For our conditions we find that, at little cost to e, 
faster detection can be obtained by omitting the effect 
of spontaneous decay in the analysis, i.e. replacing (TT]) 
by Pd — ^-iL\D{ n i)- An adaptive analysis of the same 
data set is shown in fig. [21 for t c = 500 ^s and a range 
of threshold values 10~ 6 < e c < 10 _1 . For each e c we 
plot the experimentally measured average error e against 
t a . For an average readout time of t a — 145 /is, the error 
reaches e = 1.0(1) x 10~ 4 , more than thirty times lower 
than the non-adaptive maximum likelihood method, or 
over three times as fast for the same e. The average read- 
out time for the bright (dark) state is 72 /us (219 /is). 




100 200 300 400 500 BOO 700 800 900 1000 
bin time t b [overage bin time for adaptive method A] (/is) 

FIG. 2: Average readout error e versus readout time for the 
40 Ca + optical qubit, for three different analysis methods. 
Analyses of the experimental data are plotted as symbols; 
statistical uncertainty is at most the size of the symbols. The 
accompanying curves are simulations of 10 9 trials using ideal 
Poisson statistics. Cusps arise from the discrete nature of pho- 
ton counting. T: photon-count threshold method, where the 
threshold n c was optimized for each bin time. M: maximum 
likelihood method. A: adaptive maximum likelihood method. 
Inset: Simulations using method M, showing the asymptotic 
average error e x and the time ti.i required to reach e = l.ltoo 
as functions of the net photon collection efficiency. 



For all of the analysis methods described, the bright 
S 1/2 state can be detected more accurately than the dark 
D 5/2 state because the known bright — > dark transfer rate 
is so low (we measure it to be < 10 -3 s _1 ). The data im- 
ply the bright state can be detected with over 99.9998% 
fidelity, es<2xl0~ 6 , in average time i a = 91 /is (whilst 
retaining ep, = 2.5(2) x 10~ 4 in i a — 292 /is). This asym- 
metric readout error could be exploited for QEC. 

The optical qubit offers high-fidelity readout but it suf- 
fers from two drawbacks: a finite lifetime r and the need 
for high frequency stability in the optical domain. Qubits 
stored in hyperfine ground states avoid these problems 
and exhibit some of the longest coherence times ever 
measured [5j, |8j. We turn now to the implementation 



of state detection for a qubit stored in the ||) = 



/2 

and ||) = S^/J 4 hyperfine ground states of 43 Ca + (where 
the superscripts give the quantum numbers F, Mp). For 
readout, we first map this hyperfine qubit to the optical 
qubit by state-selective transfer ||) — > D 5/2 . 

We Doppler-cool a single 43 Ca + ion on the S 3 1/2 ^> P* /2 
and Sj /2 «-> Pj /2 transitions. State preparation of \l) is 
by optical pumping with er + -polarized light on the same 
transitions. Preparation of F is achieved by extin- 
guishing the S 3 1/2 <-> P* /2 beam before the S* /2 <-» P*, 2 
beam: this prepares a statistical mixture of Mp states 
rather than |f) but the calculated effect on the measured 
readout error is < 10~ 4 . To map ||) to the D 5/2 shelf 
the ion is illuminated with 393 nm a + light which drives 
||) «-» P^J 5 . This results in optical pumping to D 5/2 
or D 3/2 with branching ratios 5.3% and 0.63% respec- 



tively 12j. Light at 850 nm containing a + and 7r polar- 
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izations is used to empty the relevant D 3/2 states via the 
P3/2 5 state (fig. [3 right). Assuming perfect repumping 
from D 3/2 , the situation reduces to that for ions without 
low-lying D levels, but with the important difference that 
we only need to drive enough transitions to transfer 1 1) to 
the shelf (rather than to collect fluorescence) before off- 
resonant excitation of |t) occurs; since the P 3/2 — » D 5/2 
branching ratio is much greater than typical photon col- 
lection efficiencies this gives a significant advantage. 

To optimize parameters the shelving transfer process 
was modelled by rate equations applied to the entire 
144-state (4S, 4P, 3D) manifold. The optimum shelv- 
ing method would be a repeated sequence of three laser 
pulses (393 nm a + , 850 nm er + , 850 nm tt), since then 
there is negligible probability of the ion decaying to 
The main limitation is off-resonant (by 3.1GHz) excita- 
tion of |t) «-> Pj)J 4 . Continuous excitation allows sim- 
ilar fidelity if the 850 nm ir component is weak (though 
with slower shelving; see fig. [3|). Accordingly, in the ex- 
periment, we used a single circularly-polarized 850 nm 
beam travelling at a small angle (2.3°) to the quantiza- 
tion axis, giving polarization components with intensities 
= (0.9992,0.0008,2 x 10~ 7 ) x 230(70)/ , 
where the saturation intensity Jo = 98 /iW/mm 2 . The 
shelving transfer was accomplished with a single simul- 
taneous 393 nm + 850 nm pulse, with duration set to the 
optimum value predicted by the model, tx — 400 fis. No 
improvement was found by varying tx or by using alter- 
nating 393 nm and 850 nm pulses (to avoid two-photon 
effects, which are not modelled by the rate equations). 
After shelving, state detection proceeds as for 40 Ca + ; 
photon-count thresholding with % = 2 ms was used as 
the optical qubit readout is not the limiting factor. 

The net average readout error for the hyperfine qubit 
was measured from 20000 trials to be e hfs = i(e-f + ejj = 
2.3(3) x 10~ 3 , a fidelity of 99.77(3)%. (The error for 
the |t) state was ej = 1.7(4) x 10 -3 .) A separate ex- 
periment measured the optical qubit readout error to be 
1.0(2) x 10 -3 , significantly poorer than for 40 Ca + due to 
a lower fluorescence rate Rb = 7500 s -1 (caused partly 
by coherent population trapping effects, which could be 
eliminated by polarization modulation techniques [13] )• 
This implies an average error ex — 1-3(4) x 10 -3 for the 
shelving transfer. This is somewhat above the modelled 
value of 0.44 x 10~ 3 ; possible reasons for this include 
imperfect circular polarization of the 850 nm beam, im- 
perfect population preparation, and broadening of the 
393 nm transition (e.g. due to finite laser linewidth). 

In conclusion, we demonstrate fast, direct, single-shot 
readout from optical and hyperfine trapped-ion qubits 
at fidelities comparable with those required for a fault- 
tolerant QC. An adaptive detection method reduced the 
optical qubit readout error by a factor m 35 for a given 
average detection time, reaching 99.99% fidelity in 145 /xs, 
with negligible time overhead due to the classical control 
system. A simple and robust method for efficient readout 




total shelving transfer time t T (sec) 



FIG. 3: Rate equation simulations for the 43 Ca + shelving 
transfer (|t), |t}) — > (|t),D 5 / 2 ). The optimum average shelv- 
ing transfer error er is shown versus total time tr allowed for 
the transfer, for both continuous and pulsed methods. The 
393 nm intensity and (in the pulsed method) all pulse dura- 
tions were numerically optimized for each tr- The error in- 
creases at short tr due to higher 393 nm intensity causing off- 
resonant shelving of |t); at long tr the error increases because 
of decay from the D 5 / 2 shelf. The minimum error is 3.5 x 10 -4 
at tr = 160 /xs; for tr = 10 /xs an error of 4.8 x 10~ 4 is avail- 
able. Right: Simplified 43 Ca + level diagram; states relevant 
to the shelving transfer are labelled by F, Mf- 

of hyperfine qubits, capable of comparable fidelity/time 
performance, was proposed and implemented. An in- 
crease in the photon collection efficiency by an order of 
magnitude, to 2%, would speed up both readout meth- 
ods, and reduce the optical qubit readout error, by a 
similar factor (fig.HJ inset). 
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